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ABSTRACT 


A  laboratory  experiment  was  conducted  examining  the  effects  of 
large  slope  surface  roughness  on  mode  propagation  in  a  rigid  walled 
waveguide.  A  waveguide  hundreds  of  depths  in  length  was  developed,  as 
well  as  a  mode  transduction  system  comprised  of  electret  type  mode 
transducers,  and  an  amplitude  and  phase  selectable  amplifier.  A  device 
for  automatically  measuring  surface  roughness  was  also  built.  The 
interaction  of  modes  with  two  types  of  surface  roughnesses, 
deterministic  and  random,  is  compared  to  propagation  in  the  smooth 
waveguide.  The  results  demonstrate  coupling  between  modes  and 
difference  of  this  coupling  in  the  presence  of  either  a 
deterministically  rough  surface  or  a  randomly  rough  surface.  A 
significant  change  in  modal  phase  speed  due  to  roughness  is 
demonstrated  at  selected  frequencies. 
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III.  MEASUREMENT  THEORY 


For  simplicity,  consider  propagation  of  a  single  frequency.  The 
method  for  measuring  modes  is  based  upon  the  observation  that  the 
pressure  at  any  point  in  the  waveguide  is  a  superposition  of  each  of 
the  modes  present  in  the  waveguide  at  that  point  [Ref.  4],  At  a  fixed 
range,  the  observed  pressure  varies  with  depth  as  the  sum  of  all  the 
modes  varies.  This  is  shown  for  an  example  of  three  modes  at  a  fixed 
instant  in  time.  (Fig.  6)  If  measurements  are  taken  at  as  many  depths 
as  there  are  modes,  a  system  of  linear  equations  is  obtained  such  as 
the  following: 


P(2q) 

=  a00 

P0 

+  a01 

P1 

+  aQ2 

OJ 

Q- 

P(z^) 

=  a1Q 

P0 

+  a1  1 

P  1 

+  a12 

P2 

P(z2) 

=  a2Q 

P0 

+  a21 

P1 

+  a22 

P2 

The  terms  on  the  left  hand  side  of  each  equation  are  the  pressures 
observed  in  the  waveguide  at  the  corresponding  depths.  The  Pn  terms  on 
the  right  hand  side  are  the  pressure  amplitudes  of  the  modes  assumed 
present  in  the  waveguide.  The  transformation  that  relates  these  two 
quantities  is  contained  in  the  coefficients  a^j,  where  i  is  the  depth 
coordinate  and  j  is  the  mode  contribution.  The  notation  employed 
relates  directly  to  the  modes.  Mode  0  is  the  plane  wave  mode  and  it 
has  the  j  subscript  0.  Mode  1  is  the  next  higher  mode,  having  pressure 
maxima  at  the  walls  but  going  through  one  minimum  at  mid  depth.  Mode  1 
has  the  j  subscript  1.  Logical  extension  shows  that  the  j  subscript, 


4.  The  Timing  Simulator 


The  Interface  Technology  model  RS-648  timing  simulator  triggers 
the  signal  generator,  OPHELEA,  and  the  wave  analyzer  at  the  intervals 
required  to  provide  coordination  between  the  system  components.  It  has 
a  rated  timing  resolution  of  50  nanoseconds. 

A  detailed  listing  of  the  equipment  settings  is  contained  in 
Appendix  C. 

Table  1 

Filtering  and  Monitoring  Equipment 

Component 
Oscilloscope 
Frequency  synthesizer 

Frequency  counter 
Bandpass  filter 

Voltage  amplifier 


Description 

Kikusui  COS5060  dual  trace 

General  Radio,  1162A  decade 
frequency  synthesizer 

Hewlett  Packard,  model  5223L 

Krohn-Hite  frequency  filter, 
model  3550 

Princeton  Applied  Research 
Amplifier,  model  113 


a.  OPHELEA 


OPHELEA  is  a  computer  designed  especially  for  data 
acquisition  and  analog/digital  conversion.  It  was  designed  by  the 
Special  Projects  Section  of  the  Naval  Air  Development  Center  with  its 
primary  component  an  Interdata  model  70  computer  that  is  programmable 
in  both  FORTRAN  and  BASIC.  It  has  a  64  thousand  byte  core  memory  and 
magnetic  disk  drives.  Data  is  stored  in  12  bit  words.  The  computer  is 
connected  to  a  Phoenix  analog  to  digital  converter,  model  ADC  712. 

This  A/D  converter  is  capable  of  encoding  input  signals  of  up  to  10 
volts  into  digital  form  with  an  accuracy  of  one  part  in  4000.  The 
system  is  extremely  fast,  with  a  maximum  data  sampling  rate  of  450,000 
samples  per  second.  Output  from  the  system  is  through  a  Texas 
Instruments  data  terminal  which  is  a  combined  keyboard  and  line 
printer. 

b.  Nicolet  FFT  Wave  Analyzer 

The  Nicolet  Scientific  Corporation  model  660B  wave  analyzer 
is  programmable  and  performs  its  own  analog  to  digital  conversion.  The 
sampling  rates  are  limited  to  manufacturer  selected  values  and  a 
maximum  of  250,000  samples  per  second.  It  has  a  built-in  CRT  screen 
that  allows  real  time  observation  of  either  the  received  signal  or  its 
spectrum. 

c.  Input  Filters,  Amplifiers,  and  Monitors 

A  summary  of  the  other  instruments  used  in  received  signal 
filtering  and  monitoring  is  included  in  Table  1. 


b.  Wavetek  Arbitrary  Waveform  Generator 

The  Wavetek  Arbitrary  Signal  Generator,  model  175,  is  a 
multifunction,  programmable  signal  generator  that  is  crystal  oscillator 
controlled.  The  crystal  oscillator  provides  excellent  frequency 
stability  and  phase  locking,  but  enables  periodic  transmission  at 
discrete  frequencies  only.  The  programmable  feature  allows  the  user  to 
drive  the  source  with  a  signal  of  arbitrary  shape.  This  is  used  to 
drive  the  source  with  a  signal  of  several  strong  Fourier  components, 
though  all  the  harmonics  are  still  limited  to  multiples  of  the  crystal 
oscillator  fundamental. 

2.  The  Receive  System 

The  receiver  system  employs  Bruel  and  Kjaer  1/2"  microphones  as 
the  primary  detection  device.  The  pre-amplifiers  are  Bruel  and  Kjaer, 
model  2619.  To  take  measurements  of  the  pressure  inside  the  waveguide, 
the  microphone  is  attached  to  a  vertical  probe  made  of  4  mm  outside 
diameter  and  3  mm  inside  diameter  steel  tubing.  This  introduces 
significant  viscous  losses  from  the  tubing  walls,  but  these  are 
minimized  by  tuning  the  probe  length  to  the  fundamental  frequency. 

Prior  to  signal  processing,  the  signal  is  amplified  by  an  HP  voltage 
amplifier  to  reduce  the  effects  of  line  noise  entering  the  system 
through  transmission  lines, 

3.  The  Signal  Processing  System 
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Figure  4.  Mode  Sources  and  Source  Frequency  Response. 


frequencies  up  to  12000  Hertz  (Hz)  and  73  percent  for  frequencies  up  to 
16000  Hz  (Fig.  4).  This  simplifies  the  analysis  of  data  for 
frequencies  that  would  otherwise  have  high  modal  densities. 

3.  Random  Surface  Measuring  Device 

The  Random  Surface  Measuring  device  (RSM)  allows  the 
accumulation  of  a  large  data  base  for  statistics  on  a  given  surface. 

It  is  capable  of  taking  1,200  samples  of  the  surface  per  hour  and  with 
a  depth  discrimination  of  20  micrometers.  This  large  sampling  provides 
meaningful  statistics  and  the  resolution  provides  a  high  degree  of 
accuracy.  Details  of  its  construction  and  operation  are  included  in 
Appendix  B. 

B.  OCEAN  PHYSICS  LABORATORY 

The  experiment  was  conducted  in  the  Ocean  Physics  Laboratory  at  the 
Naval  Postgraduate  School.  A  diagram  of  the  experimental  set-up  is 
shown  in  Figure  5  and  is  described  below: 

1.  Transmit  System 

The  mode  sources  and  multichannel  amplifier  have  already  been 
discussed  in  section  II-A-2.  For  point  sources,  Bruel  and  Kjaer 
microphones  are  used.  In  all  cases,  the  driving  signal  is  from  either 
a  Wavetek  Signal  Generator  or  a  Wavetek  Arbitrary  Waveform  Generator 
(ARB) . 

a.  Wavetek  Signal  Generator 

The  Wavetek  Phaselock  Signal  Generator,  model  186,  provides 
selectable  output  signal  forms,  such  as  a  square  wave  or  sinusoidal 
wave.  Its  primary  advantage  over  the  ARB  is  that  it  covers  a 
continuous  frequency  spectrum. 
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b.  Final  Design 

Controlled  by  these  constraints,  a  final  design  for  the 
waveguide  evolved.  The  resultant  waveguide  is  14.6  meters  long,  1.23 
meters  wide  and  0.05  meters  in  depth.  It  is  constructed  of  aluminum 
plates  and  gypsum  board.  Four  aluminum  plates,  each  3.66  meters  long 
by  1.23  meters  wide  and  0.0095  meters  thick  are  used  for  the  top 
surface  of  the  waveguide.  These  plates  are  reinforced  with  aluminum 
channel  to  minimize  their  sagging.  They  are  suspended  from  the  ceiling 
by  wires  and  turnbuckles  and  joined  together  to  make  the  14.6  meter  by 
1.23  meter  top  surface.  By  adjusting  the  turnbuckles,  the  height  of 
the  waveguide  is  maintained  at  5.19  +  .14  cm.  The  top  plates  are  mass 
loaded  by  a  layer  of  sand  to  damp  shear  waves  which  may  propagate  along 
them.  Holes  are  drilled  in  these  plates  every  0.91  meters  starting 
0.405  meters  from  one  end.  These  holes  accommodate  a  microphone  probe 
tube  of  5  mm  diameter  which  samples  the  sound  field  in  the  waveguide  at 
these  intervals.  The  bottom  surface  of  the  waveguide  is  made  of  gypsum 
board  sections  2.46  meters  long  by  1.23  meters  wide  by  0.013  meters 
thick.  The  gypsum  board  is  used  since  it  provides  a  smooth  rigid 
surface  and  does  not  support  shear  waves. 

To  achieve  the  desired  cylindrical  spreading  of  the  sound 
field,  fiberglass  absorbers  are  placed  along  the  sides  of  the  wave¬ 
guide.  Common  dense  duct  insulation  of  5  pounds  per  cubic  foot 
scatters  more  than  it  absorbs  at  shallow  grazing  angles  of  less  than  15 
degrees.  Therefore,  low  density  fiberglass  of  one  pound  per  cubic  foot 
is  used.  Signal  interference  from  side  reflections  in  the  initial 
experiments  showed  the  need  to  increase  the  distance  between  wedges 
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The  experiment  conducted  for  this  thesis  provides  some  preliminary 
data  in  this  area.  A  cylindrical  laboratory  waveguide  was  chosen  in 
that  it  provided  a  prototype  of  real  ocean  ducts,  yet  it  was  simple  and 
sufficiently  controlled  to  allow  examination  of  the  mechanism  of  mode 
interactions.  The  experiment  was  designed  to  examine  the  effects  of 
both  deterministically  .rough  surfaces  and  randomly  rough  surfaces. 


I.  INTRODUCTION 


The  purpose  of  this  thesis  is  to  examine  mode  interaction  due  to 
surface  roughness  in  a  rigid  walled  waveguide.  Of  particular  interest 
are  roughnesses  with  large  slopes  and  sharp  edges.  Though  much 
theoretical  and  experimental  work  has  been  done  for  gentle  roughnesses, 
where  the  Kirchhoff  approximation  is  applicable,  development  of  a 
theory  for  surfaces  with  large  slopes  only  began  in  the  1950's  with 
independent  work  by  Twersky  [Ref.  1]  and  Biot  [Ref.  2]  Biot  treated 
roughness  elements  as  monopole  and  dipole  scatterers  and  summed  the 
effects  of  these  scatterers  over  the  surface.  Most  notably  this  work 
has  been  continued  by  I.  Tolstoy  [Ref.  3]  who  has  extended  the  ideas  of 
Biot.  From  Tolstoy's  work  [Ref.  4]  for  example,  decreases  in  phase 
velocities  and  the  existence  of  boundary  waves  are  predicted  for 
grazing,  low  frequency  propagation  over  a  rough  surface.1  L.  Dozier 
[Ref.  5]  has  studied  the  problem  for  sonar  applications  by  using  the 
parabolic  equation  and  S.  McDaniel  [Ref.  6]  has  written  about  the 
effects  of  roughness  on  mode  coupling  in  a  waveguide.  Yet, 
developement  of  a  cohesive  understanding  of  roughness  effects  in  a 
confined  channel  has  been  hindered  by  a  lack  of  experimental  data. 


Tolstoy's  waveguide  development  is  for  a  rigid  bottom 
and  a  pressure  release  surface.  Therefore,  his  predictions 
do  not  apply  directly  to  this  thesis,  but  the  general 
phenomena  he  predicts  provide  insight  into  the  effects  being 
investigated . 
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or  mode  number,  corresponds  to  the  number  of  minimums  that  mode  has  in 


a  rigid-rigid  waveguide.  Examining  the  a^j's,  the  method  of  this 
becomes  clear.  For  example,  aoo=a10=a20=1 *  however,  aoi=“a21=1  and 
a.|i=0.  For  a  fixed  second  subscript,  j,  the  a^j's  correspond  to  the 
eigenfunction  of  the  mode  j.  Recall  that  for  a  rigid-rigid  waveguide, 
the  normalized  eigenfunctions  for  each  mode  j  are  cos(jzkz)  where  kz  is 
pi  divided  by  the  waveguide  height,  h.  Direct  comparison  of  this  form 
of  the  eigenfunction  with  the  a^j's  will  show  exact  agreement  when  z 
takes  on  values  of  0,  h/2,  and  h. 

In  order  to  obtain  the  modal  distribution  in  the  waveguide  then,  it 
suffices  to  take  measurements  of  the  pressure  at  each  depth  and  solve 
the  linear  system  P(z)=A  Pm,  where  A  is  the  matrix  of  coefficients  a^j 
(Fig.  7). 

This  simple  model  is  slightly  complicated  because  the  P(z)  column 
vector  is  complex.  It  contains  both  phase  and  amplitude  data,  but  this 
is  a  bonus.  Solving  the  linear  system  yields  amplitude  information 
about  the  mode  distribution  as  well  as  the  relative  phase  relations  of 
the  modes  present  at  the  measurement  frequency. 

The  above  description  is  most  easily  understood  for  a  single 
frequency.  In  practice,  a  non-sinusoidal  signal  is  propagated  and  the 
FFT  of  the  received  signal  is  analyzed.  Take  as  an  example  a 
triangular  pulse;  multiplying  its  frequency  spectrum  by  the  frequency 
response  of  the  electret  transducers  results  in  a  transmitted  signal 
spectrum  with  strong  components  up  to  24,000  Hz  (Fig.  8).  For  the 
smooth  waveguide,  this  one  data  run  provides  information  about  many 
frequencies  and  with  a  signal  to  noise  ratio  greater  than  20  dB.  This 
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composite  waveform  can  still  be  transmitted  with  the  primary  mode 
pressure  distribution  desired  by  using  the  mode  source.  That  is  to 
say,  propagating  a  single  mode  is  not  constrained  to  sinusoidal  signals 
only.  The  received  data  can  then  be  processed  in  the  frequency  domain. 

Processing  in  the  frequency  domain  is  done  by  the  "Modes"  progam 
(APP.  D).  The  real  and  imaginary  parts  of  the  FFT  are  transformed  by 
matrix  multiplication  with  the  inverse  of  the  matrix  A.  This  solves 
for  modal  distributions  of  both  the  real  and  imaginary  parts.  For  a 
given  mode,  the  square  root  of  the  sum  of  the  squares  of  these  two 
parts  is  the  mode’s  amplitude  and  the  inverse  tangent  of  the  ratio  of 
the  imaginary  to  the  real  part  is  the  phase. 

The  percent  energy  in  a  mode  at  any  fixed  range  is  calculated  from 
the  ratio  of  the  pressure  squared  in  that  mode  to  the  sum  of  the 
pressures  squared  in  all  modes.  This  sum  of  pressures  squared  in  all 
modes  is  proportional  to  the  total  energy  and  is  used  when  total  energy 
is  referenced. 

The  phase  information  can  be  used  to  determine  the  phase  speeds  of 
the  modes  and  any  change  of  phase  shift  due  to  roughness.  The  absolute 
phase  speeds  are  calculated  from  the  phase  change  between  measurement 
stations.  This  change,  (the  subscript  n  is  the  measurement  station 
number),  between  stations  is  taken  to  be 

^n  =  *n  *  *<n-1) 

and 

=  cuAr/Cp 

where  Ar  is  the  0.91  meters  between  measurement  positions  and  Cp  is  the 
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phase  speed.  Since  a),  A<£,  and  at  are  known,  a  value  for  Cp  can  be 
obtained.  The  phase  speed  of  mode  m  can  also  be  used  to  solve  for  the 
free  field  speed  of  sound  cQ,  using  the  equation  below. 
cQ  =  [  cp-2  +  (n/2cpfh)2  r1/2 

This  equation  is  obtained  from  the  theoretical  expression  for  the  phase 
speed  of  mode  n  in  a  rigid  walled  waveguide  of  height  h  and  frequency  f 
in  the  absence  of  absorption.  (Absorption  is  considered  later  in 
Section  V.) 

The  phase  shift  induced  by  the  roughness  can  be  calculated  from  the 
phase  speeds  for  both  the  smooth  and  rough  waveguide.  Since  the  phase 
is  given  by  the  term  <f>  =  cut-kr  and  k  =  tu/Cp,  the  roughness  induced 
phase  shift  can  be  calculated  from 

A<£  =  CJbr[(1/CpS)-(1/cpr)] 

where  the  3  subscript  denotes  the  smooth  surface  and  the  r  subscript 


the  rough  surface. 


IV.  EXPERIMENTAL  PROCEDURE 

A.  MEASUREMENTS  IN  THE  SMOOTH  SURFACED  WAVEGUIDE 

Mode  distribution  in  the  smooth  waveguide  is  the  reference  for  all 
other  measurements.  Therefore,  prior  to  making  any  meaningful 
evaluations  of  roughness  effects,  it  is  necessary  to  develop  base  line 
data  of  the  mode  distributions  in  the  smooth  waveguide.  Two  points  of 
key  interest  are  the  reproducibility  of  the  data  and  the  attenuation  of 
the  sound  field  due  to  absorption  and  divergence. 

1.  Preliminary  Validation 

A  series  of  pulsed  experiments  was  run  with  the  receiver  at 
mid-depth.  Part  of  this  experiment  was  to  confirm  the  divergence  of 
the  sound  field  observed  during  the  continuous  wave  experiments  (Fig. 
9).  Evaluation  of  the  effects  of  repeatedly  removing  the  side 
absorbers  and  then  re-inserting  them  was  also  made  during  these  pulsed 
trials.  This  is  considered  critical  in  that  the  absorbers  are 
removed  and  replaced  several  times  while  placing  the  rough  surfaces  in 
the  waveguide.  No  received  signal  variations,  other  than  those  caused 
by  background  noise,  were  observed  due  to  movement  of  the  fiberglass 
absorbers  (Fig.  10).  The  effects  of  absorber  movement  were  also 
analyzed  for. the  case  continuous  wave  transmission  and  showed  the  same 
results.  The  variations  in  the  frequency  spectrum  were  not  greater 
than  those  caused  by  background  noise,  and  this  was  less  than  0.5  dB 
re  1  volt  (Table  2). 
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Data  obtained  from  9700  Hz,  10  cycle  tone  burst 

.7  DB/M  LINEAR  LOSS  INCLUDED 


d  Attenuation  In  Smooth  Waveguide 


The  received  time  signal  of  7813  Hz  tone  burst  for  four  distinct 
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Effect  of  Absorber  Movement 


Table  2 


Effects  of  Absorber  Movement  on  Frequency  Spectrum 
for  4  Data  Runs. 


Spectrum 

level  (dB 

re  1  V) 

for  selected 

frequencies  (Hz 

Run  It 

7750 

11750 

15750 

23500 

31250 

1 

-17.4 

-16.3 

-27.0 

-36.7 

-52.0 

2 

-17.6 

-16.6 

-27.0 

-36.9 

-51  .9 

3 

-17.4 

-16.3 

-27.0 

-36.6 

-51.7 

4 

-17.2 

-16.0 

-26.5 

-36.2 

-51  .2 

Standard 

deviation  in  spectrum 

level  (dB 

re  1  V) 

7750 

11750  15750 

23500 

31250 

due  to  0.16 

movement 

(4  run  average) 

0.24  0.25 

0.29 

0.33 

due  to  0.15 

noise 

(10  run  average) 

0.23  0.21 

0.12 

0.39 

a 
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2.  Continuous  Wave  Measurements  in  the  Smooth  Waveguide 

The  continuous  wave  data  runs  consist  of  data  sets  taken  every 
0.91  meters.  A  data  set  is  comprised  of  six  or  ten  data  pairs,  which 
are  the  real  and  imaginary  parts  of  the  FFT  of  100  time  averages  of  the 
received  signal.  The  number  of  data  pairs  corresponds  to  the  number  of 
depths  at  which  the  sound  field  was  measured  at  each  range.  Whether 
six  or  ten  depth  measurements  were  used  depended  on  the  highest 
frequency  a  data  run  examined.  For  frequencies  below  16,000  Hz,  six 
depths  are  sufficient,  but  for  frequencies  up  to  31250  Hz,  ten  depths 
are  needed  because  of  the  additional  modes  which  can  propagate. 

The  signal  for  these  measurements  is  a  triangular  pulse  with  a 
fundamental  frequency  component  of  3906.6  Hz  and  strong  harmonics  up  to 
31250  Hz.  For  frequencies  below  16,000  Hz,  only  four  modes  propagate 
in  a  waveguide  of  5  centimeters  height.  Hence,  six  equally  spaced 
depth  measurements  provide  sampling  for  more  modes  than  is  required. 

For  the  smooth  waveguide,  the  energy  that  is  detected  in  these  non¬ 
propagating  modes  is  a  measure  of  the  accuracy  of  the  measurements. 
Consider  mode  5  which  is  non-propagating  for  frequencies  below  17,000 
Hz.  With  a  15,750  Hz  signal,  the  fifth  mode  is  theoretically  down  110 
dB  one  meter  from  the  source.  Hence,  it  is  reasonable  to  assume  that 
any  energy  detected  in  this  non-propagating  mode  is  due  to  measurement 
inaccuracy.  This  provides  a  qualitity  check  for  each  data  run  and 
since  there  is  typically  less  than  one  percent  energy  in  these  modes, 
the  method  is  considered  accurate  enough  to  observe  the  variations  in 
amplitude  caused  by  surface  roughnesses. 
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By  repeating  these  data  runs  several  times,  it  is  possible  to 
show  that  the  smooth  waveguide  amplitude  observations  are  quite 
reproducible.  The  standard  deviation,  over  several  days  of 
measurements,  in  the  percent  energy  present  in  a  mode  at  a  given  range 
is  typically  less  than  three  percent,  and  seven  percent  in  the  worst 
case.  This  is  demonstrated  (Table  3)  by  tabulation  of  the  energy 
distribution  in  the  smooth  waveguide  for  several  different  data  runs. 
Table  4  demonstrates  reproducibility  of  data  runs  at  the  same  location 
on  the  same  day.  At  7813  Hz  and  15625  Hz,  the  standard  deviation  in 
the  observed  percent  energy  in  any  mode  is  three  percent  or  less.  At 
11719  Hz,  the  standard  deviation  for  the  driven  mode,  mode  1,  is  seven 
percent,  but  it  is  less  than  three  percent  for  all  other  modes  at  this 
frequency.  For  frequencies  up  to  32000  Hz,  similar  arguments  can  be 
made  for  ten  depth  measurements. 

Additional  monitoring  for  erratic  or  unexpected  systematic 
response  of  the  electret  sources  is  achieved  by  using  a  1/2"  Bruel  and 
Kjaer  microphone  placed  60  centimeters  from  the  source.  The  greatest 
effect  observed  by  this  monitor  was  a  monotonic  decrease  of  six  percent 
in  the  source  amplitude  over  a  four  hour  period. 

B.  MEASUREMENTS  IN  THE  ROUGH  SURFACED  WAVEGUIDE 
1 .  Surface  Roughness  Elements 

Two  types  of  surface  roughnesses  are  used.  The  first  is  a 
deterministic,  periodic  surface  of  wedges.  The  other  surface  is 
randomly  rough,  but  of  known  statistics. 
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Table  3 

Reproducibility  of  Smooth  Waveguide 

Results . 

Energy  distributions  by 

percent  for 

four 

data 

runs,  six 

depths  and 

processing  by 

OPHELEA 

7813 

Hz  -  Smooth 

Waveguide  (CW) 

Mode  0  Energy  (%) 

Mode  1 

Energy 

(%) 

e  A  B 

C  D 

A 

B 

C 

D 

•  .  ' 

0  1 

1 

99 

97 

99 

— 

0  1 

0  1 

98 

98 

98 

98 

1  1 

1 

98 

97 

96 

- 

3  0 

0  1 

97 

98 

98 

98 

' 

1  0 

0 

98 

100 

99 

- 

6* 

0 

- 

71* 

99 

- 

10* 

3  4 

- 

82* 

97 

90 

Mode  2  Energy  (%) 

Mode  3 

Energy 

(%) 

e  A  B 

C  D 

A 

B 

C 

D 

*»  ■ 

1  2 

0 

0 

0 

0 

— 

1  2 

1  1 

0 

0 

0 

0 

•  - 

1  3 

2 

0 

0 

0 

- 

- 

0  1 

2  0 

0 

0 

0 

1 

1  0 

1 

0 

0 

0 

- 

3* 

1 

- 

6* 

0 

- 

0* 

0  1 

2* 

0 

3 

-  ■I 

< 
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Table  3  (Cont'd) 

11719  Hz  -  Smooth  Waveguide  (CW) 


Mode 

0 

Energy 

(%) 

Mode 

!  1 

Energy 

({) 

Range 

A 

B 

C 

D 

A 

B 

C 

D 

(m) 

0.9 

4 

3 

1 

- 

72 

73 

73 

- 

1  .8 

2 

1 

0 

1 

93 

93 

93 

83 

2.7 

4 

0 

1 

- 

88 

92 

92 

- 

3.7 

1 

2 

2 

2 

91 

92 

89 

62 

4.6 

7 

7 

7 

- 

85 

90 

81 

- 

5.5 

_ 

3 

4 

- 

- 

68 

94 

- 

6.4 

- 

5 

1 

4 

- 

93 

94 

89 

Mode 

2 

Energy 

(%) 

Mode  3 

Energy 

(%) 

Range 

A 

B 

C 

D 

A 

B 

C 

D 

(m) 

0.9 

1 

1 

5 

- 

26 

19 

20 

- 

1  .8 

0 

1 

2 

2 

5 

5 

4 

13 

2.7 

2 

2 

6 

- 

5 

5 

1 

- 

3.7 

1 

3 

4 

4 

7 

3 

4 

31 

4.6 

1 

0 

2 

- 

7 

3 

7 

- 

5.5 

— 

9' 

»  2 

- 

- 

6* 

0 

- 

6.4 

- 

V 

*  4 

5 

- 

1* 

1 

1 

Figure  13a.  Energy  Distribution  vs.  Range  at  15750  Hz,  Mode  1  Dr 


54 


Table  7  (Cont’d) 


Table  7 
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Energy  Distribution  for  Mode  1  Driven  System. 


Modal  energy  distributions  by  percent. 

Smooth  data  taken  at  six  depths  with  OPHELEA  processing. 
Rough  data  taken  at  ten  depths  with  Nicolet  processing. 
Signal-to-noise  ratio  greater  than  20dB  for  all  data. 


Mode  1  at  15750  Hz 


Smooth  Rubber  Wedges 


energy 

in 

mode 

(?) 

energy 

in 

mode 

(?) 

Range 

0 

1 

2 

3 

oth 

0 

1 

2 

3 

oth 

(m) 

1 .8 

1 

76 

3 

20 

0 

17 

56 

1 

20 

6 

2.1 

- 

- 

- 

- 

- 

3 

71 

6 

19 

- 

2.4 

- 

- 

- 

- 

- 

0 

61 

14 

14 

- 

6..  J 

2.7 

2 

77 

6 

12 

3 

11 

72 

3 

13 

1 

3.0 

- 

- 

- 

- 

- 

15 

76 

2 

4 

- 

*v' 

3.3 

- 

- 

- 

- 

- 

23 

57 

4 

10 

- 

3.7 

2 

77 

3 

15 

3 

27 

58 

0 

11 

4 

4.6 

3 

79 

3 

11 

4 

12 

39 

2 

33 

14 

— 

5.5 

2 

77 

4 

8 

9 

14 

49 

2 

22 

2 

6.4 

- 

- 

- 

- 

- 

4 

46 

2 

42 

6 

•  ' 

7.3 

— 

23 

39 

4 

29 

15 

Mode 

1  at 

23500  Hz 

Smooth 

Rubber 

Wedges 

energy 

in 

mode 

(?) 

energy 

in 

mode 

(?) 

Range 

0 

1 

2 

3 

oth 

0 

1 

2 

3 

oth 

*  •  ^**1 

(m) 

1 .8 

3 

70 

7 

18 

2 

1 

46 

1 

41 

11 

• - ! 

2.1 

- 

- 

- 

- 

- 

26 

5 

26 

34 

9 

"  ■  ~ 

2.4 

- 

- 

- 

- 

- 

18 

11 

41 

22 

8 

« 

2.7 

3 

65 

2 

21 

5 

14 

25 

27 

31 

3 

'  / 

3-0 

- 

- 

- 

- 

- 

14 

37 

27 

21 

1 

\\\ 

3.3 

- 

- 

- 

- 

- 

18 

48 

7 

22 

5 

3.7 

13 

62 

1 

17 

7 

13 

30 

27 

30 

0 

4.6 

9 

67 

2 

18 

4 

5 

16 

30 

31 

18 

5.5 

6 

64 

8 

16 

6 

4 

31 

16 

44 

5 

v.> 

6.4 

8 

66 

2 

18 

8 

14 

32 

6 

38 

10 

7.3 

12 

14 

33 

34 

7 

'  *  > 

52 
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V.  RESULTS  AND  CONCLUSIONS 


A.  DETERMINISTIC  SURFACE 
1.  Amplitude  Analysis 

Typical  results  for  the  mode  1  driven  system  are  tabulated  in 
Table  7  and  graphically  depicted  in  figure  13.  Several  qualitative 
observations  can  be  made  from  a  cursory  look  at  this  data.  First,  the 
average  percent  energy  in  modes,  other  than  the  driven  mode,  increases 
with  frequency  when  roughness  is  present.  This  average  is  the 
arithmetic  mean  of  the  percent  energy  over  all  ranges.  Second,  the 
variability  of  the  percent  energy  in  a  mode  is  greater  in  the  presence 
of  a  rough  surface  than  in  the  smooth  waveguide.  Third,  energy  is  not 
only  removed  from  the  driven  mode  when  roughness  is  present,  but  it  may 
reappear  in  that  mode  at  a  greater  range. 

The  average  percent  energy  of  modes  in  the  smooth  waveguide  is 
compared  to  that  in  the  wedge  roughness  waveguide  as  a  function  of 
frequency  and  kh  (Fig.  14a).  The  strength  of  interaction  is  evaluated 
by  the  percent  energy  in  modes  adjacent  to  the  driven  mode  and  is  seen 
in  this  figure  to  be  dependent  on  frequency.  To  examine  the  cause  of 
this  frequency  dependent  behaviour,  consider  the  roughness  of  the 
corrugated  rubber  surface  as  expressed  in  terms  of  the  Rayleigh  rough¬ 
ness  parameter,  kh  sinen  (Table  8). 

For  a  rigid  walled  waveguide,  the  Rayleigh  roughness  parameter  is 
nearly  constant  with  frequency.  This  is  because  the  grazing  angle  of 
the  mode,  f-»n ,  decreases  at  the  same  rate  the  wavelength  decreases  when 
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Table  6 


Mode  Cycle  Distances  (meters) 


Function 

of  Mode 

Number, 

n  and 

Frequency,  f  (Hz) 

.  f  , 
•n  1 

7750  ! 

11750  ! 

15750 

;  31250  : 

1 

.20 

.34 

.46 

.74 

2 

.06 

.14 

.22 

.46 

3 

X 

.06 

.12 

.30 

4 

X 

X 

.06 

.22 

\  ( 

all  distances  are  in  meters 
X  -  denotes  a  non-propagating  mode 


It  is  structurally  infeasible  to  place  probe  holes  this  close  on 
the  waveguide  because  of  interference  with  cross  supports.  The  alter¬ 
native  is  to  add  roughness  elements  to  the  leading  edge  of  the  rough¬ 
ness  and  move  the  source  a  distance  equal  to  the  width  of  the  added 
roughness.  This  has  the  same  affect  as  moving  the  probe  farther  from 
the  source;  the  only  change  is  the  length  of  rough  surface  the  sound 
interacts  with  before  reaching  the  receiver.  Therefore  any  changes  in 
the  phases  and  mode  distributions  is  due  to  this  greater  rough  surface 
interaction. 

Another  question  with  a  rough  surface  is  determining  the  plane  from 
which  the  waveguide  height  is  measured.  The  reference  plane  selected 
was  the  statistical  mean  height  of  the  roughness  elements.  A  priori, 
this  appears  to  be  as  good  a  selection  as  any  and  is  chosen  in  many 
texts  discussing  scattering  from  a  rough  surface. 
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Randomly  Rough  Surface  Statistics 


RMS  height  0.150  cm 

RMS  slope  80 

Correlation  length  0.148  cm 

Average  height  0.255  cm 

2.  Continuous  Wave  Measurements  in  the  Rough  Waveguide 

The  data  runs  consist  of  100  time  average  sets,  as  for  the 
smooth  waveguide,  Section  IV.2,  but  the  signal  used  is  dependent  on  the 
type  of  surface  roughness.  For  the  deterministically  rough  surface,  a 
triangular  pulse  with  a  fundemental  frequency  of  3906.6  Hz  produces  a 
signal  with  a  signal-to-noise  ratio  greater  than  20  dB  for  all  ranges. 
However,  for  the  randomly  rough  surface,  single  frequency  signals  must 
be  used  to  obtain  a  signal-to-noise  ratio  greater  than  20  dB  for  ranges 
out  to  9  meters.  The  higher  attenuation  caused  by  the  random  roughness 
scattering  both  transversely  and  along  the  waveguide  axis  makes  the 
pure  tone  signal  necessary. 

Measurements  must  be  taken  at  short  intervals  if  the  evolving 
mode  interactions  are  to  be  tracked.  One  gauge  of  the  required  inter¬ 
val  comes  from  ray  theory.  In  ray  theory,  the  cycle  distance  indicates 
how  frequently  a  ray  interacts  with  the  surface.  The  cycle  distances, 
as  predicted  by  ray  theory,  for  several  modes  at  various  frequencies 
are  given  in  Table  6.  From  the  entries  of  this  table,  a  suggested 
sampling  distance  can  be  determined.  An  interval  of  30  cm  is  selected. 
Though  this  is  greater  than  the  interval  suggested  by  Table  6,  it 
allows  a  complete  data  run  to  be  made  in  two  to  three  hours. 


lation  Functions  for  RandomLy  Rough  Surface 


a.  Deterministic  Surface 

The  deterministic  surface  is  made  of  wedge  corrugated,  hard 
rubber  matting.  The  height  of  these  wedges  is  2  mm  and  the  spacing 
between  wedges  is  3  mm  (Fig.  11).  To  put  the  roughness  to  wavelength 
ratio  in  perspective,  at  3800  Hz  kh=0.14  and  at  32150  Hz  kh=1.18;  "k" 
is  the  wave  propagation  constant  and  "h"  is  the  height  of  the  roughness 
element.  The  matting  is  placed  on  the  bottom  of  the  waveguide  in 
contiguous  one  meter  sections  aligned  so  that  the  wedges  lie  perpen¬ 
dicular  to  the  long  axis  of  the  waveguide, 
b.  Random  Surface 

This  surface  is  obtained  by  using  //2  aquarium  sand  spread 
on  the  bottom  plate  of  the  waveguide.  The  statistics  are  determined 
from  samples  made  by  the  RSM.  A  profile  of  the  surface  is  shown  in 
Figure  12a.  This  surface  data  is  processed  to  determine  the  RMS  height 
and  correlation  length.  These  two  numbers  are  then  used  to  describe 
the  surface.  The  probability  density  function  (pdf)  of  the  surface, 
which  is  obtained  from  a  histogram  of  the  surface  heights,  is  shown  in 
Figure  12b.  The  best  fit  obtained  to  this  pdf  is  indeterminate,  but 
not  Gaussian.  The  correlation  functions  for  five  of  the  six  data  runs 
taken  on  this  surface  are  shown  in  Figure  12c.  The  statistics  of  this 
surface  averaged  over  the  six  runs  is  shown  in  Table  5. 


t 
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Table  4 


Measurement  Reproducibility. 

Continuous  wave  transmission 

Six  data  runs  conducted  at  the  same  range,  on  the  same  day 
(S/N  >  20db  for  all  cases) 


7813 

Hz 

Energy 

-  in 

Mode  i 

(S) 

Trial 

0 

1 

2 

3 

4 

5 

1 

6 

84 

3 

1 

1 

0 

2 

1 

92 

6 

0 

0 

0 

3 

1 

94 

3 

1 

1 

0 

4 

1 

92 

5 

1 

1 

0 

5 

1 

91 

7 

0 

0 

0 

6 

3 

91 

5 

0 

0 

0 

17719 

i  Hz 

Energy  in 

Mode 

(X) 

Trial 

0 

1 

2 

3 

4 

5 

1 

11 

64 

8 

15 

2 

0 

2 

3 

81 

2 

12  - 

0 

1 

3 

5 

79 

3 

12 

1 

0 

4 

6 

83 

0 

8 

2 

1 

6 

7 

30 

0 

12 

1 

0 

15625  Hz 

Energy  in 

Mode 

(%) 

Trial 

0 

1 

2 

3 

4 

5 

1 

13 

54 

2 

19 

11 

1 

2 

10 

54 

1 

26 

9 

0 

3 

10 

53 

4 

28 

6 

0 

4 

7 

53 

2 

27 

10 

1 

Table  3(Cont'd) 


15625  Hz  -  Smooth  Waveguide  (CW) 


Mode 

0 

Energy 

(?) 

Mode  1 

Energy 

(2) 

Range 

A 

B 

C 

D 

A 

B 

C 

D 

(m) 

0.9 

2 

0 

1 

- 

66 

73 

67 

- 

1 .8 

1 

1 

2 

- 

72 

77 

78 

- 

2.7 

7 

0 

1 

- 

72 

79 

80 

- 

3.7 

6 

4 

1 

- 

76 

78 

78 

- 

4.6 

8 

5 

0 

- 

74 

81 

84 

- 

5.5 

- 

12 

1 

- 

- 

71 

32 

- 

6.4 

- 

9 

2 

- 

- 

42 

49 

- 

8.2 

- 

2 

9 

- 

- 

82 

73 

- 

Mode 

2 

Energy 

(?) 

Mode  3 

Energy 

(X) 

Range 

A 

B 

C 

D 

A 

B 

C 

D 

(m) 

0.9 

1 

0 

9 

- 

30 

25 

22 

- 

1 .8 

2 

5 

0 

- 

23 

16 

18 

- 

2.7 

9 

4 

7 

- 

10 

17 

10 

- 

3.7 

0 

0 

9 

- 

18 

17 

10 

- 

4.6 

1 

0 

9 

- 

14 

13 

7 

- 

5.5 

- 

2 

6 

- 

— 

4 

10 

— 

6.4 

- 

1 

4 

- 

- 

43 

42 

- 

8.2 

0 

9 

11 

7 

Date  of  run: 

A  - 

9/17  B 

-  9/29 

C  -  9/30 

D  - 

■  10/28 

*  -  S/N  <  20  dB 


Percent  Energy  in  Mode 


The  average  percent  energy  is  the  arithmetic  mean  over  all 
ranges  for  that  frequency. 

oo  Mode  0 

Cl 


£  Frequency  (kHz) 


c  Mo  de  2 


8  16  24  32 


Frequency  (kHz) 

Mode  3 


8  16  24  32 


Frequency  (kHz) 

0  -  smooth  surface  A  -  wedge  roughness 

Figure  14a.  Average  Percent  Energy  in  Modes  Adjacent  to  Mode  1 
(Mode  1  Driven) . 


frequency  increases.  Therefore,  the  Rayleigh  roughness  parameter  is 
not  helpful 

Table  3 

Rayleigh  Roughness  Parameter,  R,  h  =  2  mm 

Mode 

12  3  4 

R  0.13  0.25  0.38  0.55 

in  accounting  for  the  correlation  between  increasing  frequency  and 

increasing  internal  energy  transfer  due  to  the  surface  roughness.  This 

is  expected  because  the  Rayleigh  roughness  parameter  was  developed  for 

the  special  case  of  small  slopes  and  specular  scatter.  However,  as 

noted  in  section  IV,  0.12  <  kh  <  1.10  is  in  the  region  of  Rayleigh 

scatter  within  which  scattered  power  increases  as  the  fourth  power  of 

the  frequency.  It  is  noted  that  the  variation  of  energy  in  modes 

adjacent  to  the  driven  mode,  as  measured  by  the  standard  deviation 

about  the  average  percent  energy  for  all  ranges  at  a  given  frequency, 

is  also  a  strong  function  of  frequency  (Fig.  14b). 

An  interesting  effect  is  the  apparent  oscillation  of  energy 
from  one  mode  to  another.  At  31250  Hz.  for  example,  energy  is  removed 
from  mode  1  and  distributed  to  closely  adjacent  modes  as  compared  to 
the  energy  distribution  in  the  smooth  waveguide  at  2  meters  (Fig.  13c). 
By  3  meters,  the  energy  distribution  has  returned  to  one  more  like  the 
smooth  waveguide.  Oscillations  like  this  continue  throughout  the 
ranges  of  observation. 

It  is  noted  that  the  attenuation  rate  for  the  corrugated  sur¬ 
face  is  slighty  less  than  for  the  smooth  surface  (Fig.  15).  A 
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Figure  15a. 


Total  Energy  Attenuation  at  15750  Hz,  Mode  1  Source, 
Smooth  vs.  Deterministic  Roughness. 
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Figure  15b.  Total  Energy  Attenuation  at  23500  Hz,  Mode  1  Source, 
Smooth  vs.  Deterministic  Roughness. 
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possible  explanation  may  be  that  the  viscous  losses  caused  by  the 
gypsum  board  surface  may  be  greater  than  the  effects  of  the  rubber 
wedges.  A  future  experiment  to  test  this  would  be  to  use  smooth  rubber 
of  the  same  material  as  the  wedges  for  the  bottom  surface  in  the 
smooth  waveguide  reference  runs.  Regardless,  the  corrugated  wedges 
can  be  described  as  a  one  dimensional  roughness  and  thus  there  is 
minimal  scattering  of  sound  to  the  sides  of  the  waveguide.  Thus  a 
large  change  in  the  e  ergy  attenuation  rate  is  not  expected  and  is  not 
observed.  The  remaining  energy  attenuation  (dB/m),  after  subtracting 
out  attenuation  due  to  geometrical  speading,  for  mode  1  driven  in  both 
the  rough  and  smooth  walled  waveguide  is  giveri  in  Table  9.  The  energy 
attenuation  for  a  given  run  is  the  average  of  the  incremental 
attenuations,  AI, 


AI  =  -10  log 

/  Total  Energy  \ 

/first  range\ 

Distance 

iTotal  Energy  first  range] 

\  range  ) 

where  distance  =  range  -  first  range.  In  table  9,  the  values  are  the 
averages  over  three  runs. 

Similar  results  are  observed  for  the  mode  2  driven 
system  as  are  observed  for  the  mode  1  driven  system  (Table  10  and 
Figure  16).  The  average  energy  in  modes  immediately  adjacent  to  mode  2 
and  their  variation  in  energy,  as  measured  by  the  standard  deviation  of 
the  energy  about  its  mean,  increase  with  frequency  (Fig.  17).  This 
comparison  reveals  no  clear  information  about  mode  0  though. 

Oscillations  of  the  energy  in  modes  is  also  observed.  To 
determine  if  this  is  a  real  phenomenon  or  due  to  different  attenuation 


rates  for  each  mode,  the  modal  amplitudes  are  compared  as  a  function  of 
range  (Fig.  18,  19,  20).  In  these  figures,  the  rough  waveguide 


Table  9 


Energy  Attenuation 

in  Smooth  and  Wedge  Roughness  Waveguide, 
Mode  1  Driven  (dB/m). 


Freq  (Hz) 

7813 

11719 

15626' 

19531 

23438 

27344 

31250 

Smooth 

Waveguide 

1.66 

2.21 

2.67 

3.06 

3.27 

3.88 

4.15 

Wedge 

Roughness 

1.37 

2.11 

2.27 

2.36 

3.43 

2.85 

2.60 

Table  10 


Energy  Distribution  for  Mode  2  Driven  System, 
Smooth  and  Deterministically  Rough  Surfaces. 


Modal  energy  distribution  by  percent  for  ten  depth 
measurements  taken  by  the  Nicolet  wave  analyzer 
and  processed  on  the  IBM  3030. 
Signal-to-noise  ratio  greater  than  20  dB  for  all  data. 


Mode  2  at  8000  Hz 

Smooth  Rubber  Wedges 

energy  in  mode  (})  energy  in  mode  (%) 


Range 

(m) 

0 

1 

2 

3 

0 

1 

2 

3 

1  .8 

23 

1 

76 

0 

17 

1 

81 

0 

2.7 

24 

4 

72 

0 

11 

4 

80 

0 

3.7 

11 

3 

85 

0 

16 

3 

79 

0 

4.6 

22 

20 

57 

0 

19 

6 

74 

0 

5.5 

31 

13 

54 

0 

10 

6 

74 

1 

6.4 

40 

9 

51 

0 

40 

18 

38 

1 

Mode  2  at 

Smooth 

energy  in  mode  {%) 

13000  Hz 

energy 

Rubber  Wedges 
in  mode  (.%) 

Range 

0 

1 

2 

3 

0 

1  2 

3 

(m) 

1.8 

20 

3 

74 

4 

15 

2  79 

2 

2.7 

23 

5 

71 

1 

37 

2  42 

18 

3.7 

24 

1 

75 

0 

21 

1  77 

1 

4.6 

27 

4 

68 

0 

46 

1  43 

8 

5.5 

23 

2 

75 

0 

18 

7  66 

3 

6.4 

31 

3 

64 

1 

7 

2  73 

18 

7.3 

15 

4 

80 

1 

1 

4  95 

0 

8.2 

20 

3 

75 

1 

7 

12  78 

1 

H* 


65 


Table  10  (cont'd) 


Mode  2  at  15750  Hz 


Smooth  Rubber  Wedges 


energy 

in  mode 

(X) 

energy 

in  mode 

(X) 

Range 

(m) 

0 

1 

2 

3 

0 

1 

2 

3 

1.8 

13 

2 

83 

2 

5 

11 

76 

6 

2.7 

11 

1 

84 

0 

14 

12 

68 

2 

3.7 

11 

6 

79 

1 

17  • 

19 

58 

6 

4.6 

12 

7 

80 

1 

11 

13 

61 

12 

5.5 

12 

3 

83 

1 

6 

19 

73 

1 

6.4 

16 

1 

81 

0 

1 

23 

74 

0 

Mode  2  at  23500  Hz 


Smooth  Rubber  Wedges 


energy 

in  mode 

(X) 

energy 

in  mode 

(X) 

Range 

(m) 

0 

1 

2 

3 

0 

1 

2 

3 

1.8 

7 

4 

79 

5 

4 

6 

80 

4 

2.7 

11 

6 

77 

3 

5 

5 

66 

13 

3.7 

11 

19 

67 

2 

6 

14 

75 

0 

4.6 

18 

4 

74 

1 

11 

7 

68 

12 

5.5 

14 

6 

74 

4 

5 

3 

38 

1 

6.4 

20 

2 

70 

6 

6 

18 

65 

4 

Mode  2  at  31250  Hz 


Smooth  Rubber  Wedges 


energy 

in  i 

node 

(X) 

energy 

in  mode 

(X) 

Range 

(m) 

0 

1  2 

3 

4 

oth 

0 

1  2 

3 

4 

oth 

1.8 

11 

1  72 

5 

6 

5 

8 

32  17 

3  21 

19 

2.7 

11 

8  72 

2 

4 

3 

4 

13  52 

7  14 

10 

3.7 

3 

5  68 

12 

5 

7 

2 

30  34 

2  17 

15 

4.6 

3 

6  76 

8 

4 

5 

4 

34  22 

6  12 

22 

5.5 

8 

17  66 

1 

3 

1 

11 

3  49 

14  12 

11 

6.4 

22 

4  55 

7 

4 

8 

15 

3  20 

24 

8 

30 

66 
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Figure  21a.  Total  Energy  Attenuation  at  15750  Hz,  Mode  2  Source, 
Smooth  vs.  Deterministic  Roughness. 
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amplitudes  have  been  normalized.  The  normalization  factor  is  the 
square  root  of  the  ratio  of  total  energy  energy  present  in  the  smooth 
waveguide  to  total  energy  present  in  the  rough  waveguide  at  that  range. 
This  compensates  for  any  amplitude  variations,  as  much  as  2  dB,  caused 
by  the  direction  the  probe  tube  faces.  In  both  the  smooth  and  rough 
walled  waveguide,  the  rate  of  attenuation  for  different  modes  is  not 
the  same.  As  this  does  not  affect  the  distribution  of  energy  among 
modes  in  the  smooth  waveguide,  it  does  not  explain  oscillations 
observed  in  the  rough  waveguide.  The  data  suggests  that  there  is 
oscillatory  mode  coupling  in  the  presence  of  a  deterministically  rough 
surface. 

The  rate  of  total  energy  attenuation  in  the  mode  2  driven 
system  for  the  deterministically  rough  waveguide  (Fig.  21)  is  less  than 
for  the  smooth  waveguide,  as  in  the  case,  of  the  mode  1  driven  system. 
This  is  seen  in  Table  11  which  lists  the  same  information  as  Table  9 
except  mode  2  is  driven.  The  values  listed  are  averaged  over  2  runs 
for  the  rough  data  and  one  run  for  the  smooth  data.  It  is  noted  that 
these  attenuation  rates  are  slightly  different  from  those  listed  in 
Table  9  for  mode  one  driven.  Since  the  driven  mode  contains  most  of 
the  energy  present,  this  difference  can  be  interpreted  as  a  difference 
in  attenuation  for  modes  one  and  two.  From  a  ray  point  of  view  and 
assuming  the  same  loss  per  bounce,  this  is  to  be  expected  because  mode 
two  has  a  shorter  skip  distance  and  therefore  interacts  with  the 
surface  more  often  per  unit  distance.  The  difference  in  energy  levels 
observed  on  these  two  sets  of  figures  is  because  a  different  scaling 


MODE  2  SOURCE  -  WEDGE  ROUGHNESS 


RANGE  (M) 


Figure  20c.  Mode  0  Amplitude  at  31250  Hz,  Mode  2  Source, 
Smooth  vs.  Deterministic  Roughness. 
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Figure  20b.  Mode  1  Amplitude  at  31250  Hz,  Mode  2  Source, 
Smooth  vs.  Deterministic  Roughness. 
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Figure  20a.  Mode  2  Amplitude  at  31250  Hz,  Mode  2  Sourc 
Smooth  vs.  Deterministic  Roughness. 
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Figure  19c.  Mode  0  Amplitude  at  23500  Hz,  Mode  2  Source 
Smooth  vs.  Deterministic  Roughness. 
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Figure  19a.  Mode  2  Amplitude  at  23500  Hz,  Mode  2  Source, 
Smooth  vs.  Deterministic  Roughness. 
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Figure  18c.  Mode  0  Amplitude  at  15750  Hz,  Mode  2  Source 
Smooth  vs.  Deterministic  Roughness. 
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Figure  21b.  Total  Energy  Attenuation  at  23500  Hz,  Mode  2  Sovrc 
Smooth  vs.  Deterministic  Roughness. 
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MODE  2  SOURCE  -  WEDGE  ROUGHNESS 


Figure  21c.  Total  Energy  Attenuation  at  31250  Hz,  Mode  2  Source 
Smooth  vs.  Deterministic  Roughness. 


Table  11 


Ener"-,  Attenuation 

in  Smooth  and  Wedge  Roughness  Waveguide 
Mode  2  Driven  (dB/m). 


Freq  (Hz) 

7750 

15750 

23500 

31 

Smooth 

Waveguide 

2.28 

2.29 

2.69 

3 

Wedge 

Roughness 

2.11 

1 .22 

2.24 

3 

factor  was  used  for  OPHELEA  (Fig.  15  data)  than  was  used  for  the 
Nicolet  (Fig.  21  data). 


2.  Phase  Analysis 
a.  Mode  1 

The  measured  phase  of  the  driven  mode  can  be  used  to 
determine  the  phase  speeds  of  the  modes  and  any  relative  shift  in  them 
due  to  the  introduction  of  the  roughness.  Only  phase  information  for 
the  driven  mode  is  analyzed  since  that  moda  has  the  largest  amplitude 
and  so  provided  the  most  reliable  data.  For  a  given  frequency  the 
phase  speed  for  a  mode,  m,  is  given  by 

cp  =  c0/[1-(c0m/2hf)2]1/2 

where  f  is  the  frequency  in  hertz,  cQ  is  the  free  field  speed  of  sound 
in  meters  per  second,  m  is  the  mode  number  and  h  is  the  height  of  the 
waveguide  in  meters.  The  most  obvious  difference  between  the  phase 
speeds  in  the  smooth  versus  rough  walled  waveguide  is  due  to  a  change 
in  the  mean  height  of  the  waveguide.  The  mean  height  of  the  smooth 
waveguide  (0.053  m)  is  greater  than  the  mean  height  of  the  rough  walled 
waveguide  (0.051  m),  the  difference  being  the  mean  height  of  the  rough¬ 
ness  (0.051  m).  Since  the  phase  speed  is  inversely  proportional  to  the 
mean  height,  the  expected  phase  speed  in  the  rough  walled  waveguide  is 
greater  than  in  the  smooth  waveguide.  Using  these  heights  and  a  free 
field  speed  of  sound  of  3^5  m/s,  the  phase  speeds  in  the  smooth,  cn„, 

ps 

and  rough  walled  waveguide,  cpr,  are  calculated  (Table  12).  Based  on 
the  relation  =  cut-kr,  where  k  =  u/ cp,  the  phase  change  between 
measurement  positions  in  the  smooth  walled  waveguide  should  be  greater. 
The  phase  gain  (relative  spatial  phase  shift)  of  the  smooth  over  the 


7813  Hz 


11719  Hz 


15625  Hz 


Cpg  (m/s) 

379.50 

359.13 

352.74 

cpr  (m/s) 

382.72 

360.33 

353.38 

A<t>h  (°) 

56.86 

35.68 

26.34 

Table  12.  Effect  of  Waveguide  Height  on  Phase  Speed. 

A4>h  for  a  Travel  Distance  of  0.91  m  (Theoretical). 


rough  waveguide  (Table  12)  is  calculated  from 
^<f>h  =  uur[(1/cps)-(1/cpr)] 

where  ^r  is  the  distance  between  measurement  positions  (0.91m). 

is  nearly  invariant  with  respect  to  slight  changes  in 
the  free  field  speed  of  sound,  with  the  greatest  changes  occurring  at 
frequencies  close  to  cutoff.  The  stability  of  with  changes  in 
sound  speed  is  calculated  for  a  narrow  range  of  speeds  at  7813  Hz 
(Table  13). 

Another  factor  affecting  phase  speed  is  attenuation.  For 
mode  0  in  a  duct  in  which  viscous  losses  dominate  the  total  attenua¬ 
tion,  Pierce  [Ref.  8]  gives  the  apparent  phase  speed  as 

vph  =c-  (c^a/W) 

Using  a  s  a/8.7  nepers/meter,  where  a  is  the  attenuation  in  dB/meter 
given  in  Section  V-A.1,  vph  is  calculated  for  0=345  m/s  (Table  14). 

While  the  above  formula  is  for  mode  0  rather  than  mode  1  in 
a  duct  different  from  the  waveguide  used  in  this  experiment,  it  does 
give  an  estimation  of  the  relative  importance  of  dispersion  due  to 
attenuation.  Because  the  attenuation  is  nearly  the  same  in  the  smooth 
waveguide  as  it  is  in  the  waveguide  with  corrugated  roughness,  disper¬ 
sion  due  to  attenuation  can  account  for  very  little  of  the  relative 
phase  shift  between  the  smooth  and  rough  waveguide.  The  maximum 
difference  in  relative  phase  speed  between  the  smooth  and  corrugated 
rough  surface  waveguide  for  mode  one  due  soley  to  attenuation  is  0.02%. 
The  differences  in  phase  speeds  due  to  changes  in  height  therefore 
dominate  over  those  due  to  attenuation,  although  even  this  cannot 
explain  the  changes  in  the  dispersion  due  to  the  roughness. 
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c  (m/s) 

342 

343 

344 

345 

346 

Cps  (m/ s ) 

375.51 

376.84 

378.17 

379.50 

380.83 

cpr  (m/s) 

378.64 

380.00 

381.36 

382.72 

384.09 

A4»h  (°) 

56.46 

56.60 

56.73 

56.86 

57.17 

Relative  Phase  Shif 
eoretical) . 


— 

7813  Hz 

11719  Hz 

15625  Hz 

vph  in 

Smooth 

Waveguide 

344.54 

344.60 

344.62 

vph  in 

Rough 

Waveguide 

344.61 

344.61 

344.68 

The  observed  values  of  the  phase  shifts,  A0obn  have  been 
calculated  from  the  experimental  data.  If  </>rnis  defined  as  the  phase 
measured  in  the  rough  waveguide  at  position  n  (n  x  0.91m  from  the 
source)  and  <£gnis  defined  as  the  phase  measured  in  the  smooth  waveguide 
at  position  n,  the  measured  A<£obn  over  one  range  increment  is  then 

A^obn  =  ^rn“^sn^  "  ^r(n-1  )-<^s(n-1  )^ 

These  values  are  calculated  for  four  different  data  runs  for  the 

frequencies  7813  Hz,  11719  Hz,  and  15625  Hz  (Tables  15,  16,  and  17). 

The  average  (O<£obn>)  and  standard  deviation  for  each  frequency  are 
computed  and  are  also  listed  in  these  tables.  These  standard 
deviations  are  within  those  expected  for  a  signal  to  noise  ratio  of  20 
dB,  the  minimum  accepted  for  reliable  data.  With  a  20  dB  signal  to 
noise  ratio,  the  maximum  phase  error  (when  signal  and  noise  phases  are 
in  quadrature)  is  tan-1 1/10  =  5.7  degrees.  The  difference  of  two 
measured  phases  has  a  possible  error  of  11.4  degrees,  but  A<fiobn  is  the 
difference  of  two  differences  and  so  the  possible  error  induced  from 
noise  alone  is  22.8  degrees.  Data  points  which  are  in  obvious 
disagreement  with  the  others  at  the  same  frequency  are  not  included  in 
the  calculation  of  the  average  and  are  in  parentheses  in  the  tables. 
Data  points  with  an  asterisk  are  1/2  the  phase  shift  calculated  over 
1.82  meters.  This  is  necessary  in  some  instances  because  of  a  missing 
reading  at  a  particular  range. 

Upon  comparison  of  ±4>obn  with  A<£b  for  each  frequency,  it  is 
seen  that  in  all  cases  they  differ  by  more  than  the  standard  deviation 
of  If  the  assumption  is  made  that  the  roughness  causes  a 

decrease  in  phase  speed,  the  minimum  change  in  phase  speed  necessary  to 


17  SEP 


29  SEP 


30  SEP 


28  OCT 


A<frob3 

12 

31 

30 

16 

A4>ob4 

42 

-22 

44 

-1 

A<t>ob5 

18 

32 

(289) 

50 

A<t>ob6 

47 

51 

21 

38 

A*ob7 

41 

12 

-28 

A<f>ob8 

- 

- 

35 

A(^ob9 

9* 

17* 

29 

17  SEP 


29  SEP 


30  SEP 


28  OCT 


A<t>ob3 


A<>ob4 


A*ob5 


A<fob6 


A<f*ob7 


A<t>ob8 


A<J>ob9 


A(t>oblO 


-45 

-72 

-37 

-65 

-49 

-15 

-34 

-39 

-19 

-32 

-72 

-39 

-12 

-65 

-45 

-27* 

-33* 

(-119) 

-79 

<A*obn>  “ 

a  -  18.30° 


Table  16.  Observed  Relati' 
Mode  1  at  11719 


9 


<Mobn>  = 

a  =  22.62‘ 


Table  17.  Observed  Relative  Phase  Shifts  in  Degree  for 
Mode  1  at  15625  Hz. 


account  for  the  observed  phase  shifts  (<-^0bn^  can  be  computed.  If  a 
free  field  speed  of  sound  of  345  m/s  is  assumed,  cps  from  Table  12  can 
be  used  for  each  frequency.  The  observed  phase  speed  in  the  rough 
waveguide,  cpr.f  can  then  be  calculated  from 

<^obn>  =  0Mr[(1/cps)  -  (1/cpr,)] 

The  results  of. this  calculation  are  listed  in  Table  18.  The  difference 
between  the  expected  and  observed  phase  speed  in  the  rough  walled 
waveguide  can  then  be  computed  as  a  percentage  of  the  expected  value. 
(Table  18)  It  is  apparent  that  the  roughness  causes  a  change  in  the 
phase  speed  of  mode  1  that  can  not  be  accounted  for  by  the  change  in 
waveguide  height. 

As  a  check  of  the  percent  differences  obtained  above,  the 
phase  data  is  processed  by  a  second  method.  The  absolute  phase  speeds 
are  calculated  from  the  average  phase  change  between  measurement 
stations.  The  phase  change  between  stations  is  taken  to  be 

sn  =  ^sn  “  ^s(n-l) 

and 

A<Vn  "  ^rn  “  ^r(n-l) 

for  the  smooth  and  rough  waveguide,  respectively.  For  these 
calculations  only  the  data  of  September  29  and  30  are  used  since  it 
shows  the  most  consistency.  Also,  the  data  for  the  rough  waveguide  at 
15625  Hz  showed  no  consistency  and  is  not  used.  The  values  of  $<t>sn  and 
■\4>rn  are  listed  in  Tables  19,  and  20.  Their  averages  (<A«£sn>  and  <A<£rn>) 
and  standard  deviations  are  listed  in  Table  21.  This  measured  phase 
shift  is  equal  to 

A<#>  =  UMr/cp 
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UNCLASSIFIED 
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NL 


Ai*,s6 


a<*>s7 


a4>s9 


A$slO 
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29  SEP 

7813  Hz 

30  SEP 

1 

-6 

1 

15 

18 

11 

2 

3 

8 
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15625  Hz 
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84 
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61 

79 

0 
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118 
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30  SEP 
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11719  Hz 
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33 
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6 
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Since  u>,  d</>,  and  dr  are  known,  a  value  for  Cp  can  be  obtained  and  the 
phase  speeds  calculated  for  the  smooth  and  rough  waveguide  (cpsob  and 
Cprob^  (Table  22).  The  phase  speed  in  the  smooth  waveguide  is  then 
used  in  the  previously  stated  equation  for  phase  speed  to  obtain  the 
free  field  speed  of  sound,  cQ  (Table  22).  This  value  for  cQ  is  used  to 
calculate  the  theoretical  phase  speed  in  a  waveguide  of  0.051  m  height 
(cpr,,  also  in  Table  22).  The  percent  difference  between  cprob  and 
cpr,  is  calculated  (Table  22).  When  the  percent  difference  in  Table  22 
is  compared  to  the  percent  difference  listed  in  Table  18,  it  is  seen 
that  they  are  in  agreement.  This  comparison  demonstrates  that  the  data 
supports  a  change  in  phase  speed  due  solely  to  the  roughness. 

It  is  interesting  to  note  that  the  calculated  free  field 
sound  speed,  c0,  (Table  22)  increases  with  frequency.  Since  cQ  is  in 
reality  the  phase  speed  of  mode  zero,  an  increase  with  frequency  is 
what  is  expected  in  an  attenuating  medium.  The  increase  observed, 
however,  is  greater  that  that  expected  for  the  measured  attenuation. 

There  are  some  qualifications  on  these  conclusions.  First 
is  that  all  of  the  data  used  in  the  above  analysis  was  collected  with 
the  use  of  OPHELEA.  It  i3  noted  that  in  data  collected  by  OPHELEA,  the 
total  energy  of  the  modes  measured  repeatedly  at  any  one  position  under 
identical  conditions  decreased  with  time.  It  is  because  of  this  that 
OPHELEA's  reliability  was  brought  into  question  and  its  use  suspended. 
This  drift,  however,  does  not  seem  to  affect  the  measured  phase  of  the 
modes.  Test  runs  at  one  position  separated  by  two  hours  (the  approxi¬ 
mate  time  of  a  data  run)  showed  a  measured  phase  shift  of  less  than  8 
degrees.  This  8  degree  shift  is  only  about  two  degrees  more  than  that 


which  can  be  explained  by  noise  as  discussed  earlier.  It  is  felt  that 
any  phase  shift  introduced  by  the  equipment  does  not  significantly 
affect  the  above  results  since  the  results  are  consistent  from  run  to 
run.  This  assumption  remains  to  be  verified  however.  This  can  be  done 
by  taking  data  in  the  smooth  waveguide,  then  the  rough,  and  then  the 
smooth  again.  If  the  apparent  change  in  phase  speed  is  real,  compari¬ 
son  of  each  smooth  run  to  the  rough  should  yield  the  same  results. 

A  second  qualification  to  the  above  analysis  is  that  the 
assumption  is  made  that  the  phase  speed  in  the  rough  waveguide  is 
slowed  by  the  minimum  amount  necessary  to  account  for  the  observed 
phase.  If  the  number  of  complete  cycles  the  mode  underwent  is 
different  from  assumed,  the  phase  speed  of  the  mode  may  decrease  more 
than  listed  or  increase.  Tables  23,  24,  and  25  show  the  phase  speeds 
and  percent  change  in  phase  speeds  obtained  from  the  first  method 
discussed  if  360  q  degrees  is  added  to  the  observed  phase  change 
(0</>obn>)  for  various  values  of  q.  The  tables  are  only  calculated  to 
the  first  value  of  the  percent  change  greater  than  five  percent, 
similar  results  can  be  obtained  using  the  second  method  if  the  assumed 
number  of  cycles  that  mode  one  in  the  rough  walled  waveguide  undergoes 
between  stations  is  changed  by  q.  The  actual  amount  of  phase  change 
can  be  determined  if  measurements  of  the  sound  field  are  taken  at  much 
shorter  intervals, 
b.  Mode  2 

The  data  for  mode  2  is  not  as  consistent  as  that  for  mode 
1.  The  only  frequency  for  which  the  standard  deviation  of  is 

les3  than  30  degrees  is  23500  Hz.  -^o^n  at  23500  Hz  for  the  two  runs 
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N  V 


c 


c 


cprf  (m/s) 


jr  pr ' 


x  100(%) 


-2 

334.53 

5.34 

-1 

342.55 

3.07 

0 

350.96 

0.69 

1 

359.80 

-1.82 

2 

369.09 

-4.45 

3 

378.88 

-7.22 

se  Speed  of  Mode  1  over  wedge 
15625  Hz  Due  to  Cycle  Error  q. 


IH? 


V  ^  -V  V 


made  with  the  mode  2  source  and  the  wedge  roughness  are  listed  in  Table 
26  along  with  O<£obn>  and  the  standard  deviation.  The  first  method  of 
analysis  performed  on  the  mode  1  data  is  performed  on  this  data.  The 
results  for  cQ  assumed  as  3^5  m/s  are  listed  in  Table  27.  As 
previously  stated,  the  percent  change  listed  would  be  different  if  the 
actual  relative  phase  shift  are  360  q  degrees  different  than  that 
measured. 

B.  RANDOMLY  ROUGH  SURFACE 
1 .  Amplitude  Effects 

The  most  notable  effect  of  the  randomly  rough  surface  is  the 
increased  attenuation  rate  of  sound  amplitude.  The  oscillatory  effects 
observed  in  the  deterministically  rough  surface  are  not  strongly 
observed  and  there  is  a  much  clearer  coupling  of  energy  from  higher 
modes  to  lower  modes. 

The  dramatic  increase  in  energy  attenuation  caused  by  the 
randomly  rough  surface  is  illustrated  by  table  28  where  the  attenuation 
rate  AI  in  dB/m,  excluding  geometrical  spreading,  is  determined  as 
follows: 


AI  =  -10  log 

/  Total  Energy  \ 

/first  rangeN 

Distance 

vTotal  Energy  first  range) 

\  range  J 

where  distance  =  range  -  first  range.  Recall  that  for  the 
deterministically  rough  surface,  the  energy  attenuation  rates  in  the 
smooth  and  rough  waveguide  were  nearly  identical  (Tables  9,  11).  This 
same  comparison  is  presented  for  the  randomly  rough  surface  with  mode  1 
driven  (Fig.  22)  and  accents  the  greatly  increased  attenuation  caused 


by  the  randomly  rough  surface. 


7750  Hz 

15750  Hz 

31250  Hz 

315 

194 

251 

(165) 

A<?ob4 

302 

211 

174 

257 

Al?ob5 

313 

203 

274 

227 

A*ob6 

300 

209 

220 

252 

A<?)ob7 

310 

191 

236 

274 

A<fob8 

313 

220 

242 

244 

A4>ob9 

314 

210 

254 

^A^obn^ 

310 

205 

'229 

231 

0 

6.05 

10.12 

32.56 

■■ 

Table  34.  A4>obn  and  <A<J>obn>  in  Degrees  for  Mode  1  Driven 
over  Randomly  Rough  Surface. 
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7750  Hz 

15750  Hz 

23500  Hz 

31250  Hz 

cps  (m/s) 

380.15 

352.61 

348.36 

346.89 

cpr  (m/s) 

384.44 

353.43 

348.71 

347.08 

A$h  (°) 

74.69 

34.02 

22.23 

16.19 

Table  33.  Phase  Speeds  in  Che  Smooth  and  Randomly  Rough  Waveguide  for 
Mode  i  Based  upon  a  Free  Field  Sound  Speed  of  345  m/s. 

A<J>h  for  a  Travel  Distance  of  0.91  m  (Theoretical). 
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rough  surface  is  changing  the  boundary  condition  toward  pressure 
release.  It  is  not  expected  that  modes  0  and  1  would  remain  out  of 
phase  at  the  bottom  of  the  waveguide  at  several  consecutive  ranges  for 
two  frequencies,  so  an  altered  boundary  condition  due  to  surface 
roughness  is  possible. 

2.  Phase  Analysis 

Due  to  the  fact  that  the  randomly  rough  surface  could  not 

easily  be  laid  down  or  removed,  there  is  only  one  data  set  for  which 

runs  with  the  randomly  rough  walled  waveguide  and  smooth  waveguide  were 

made  on  the  same  day.  Phase  information  was  analyzed  by  the  first 

method  discribed  in  Section  V-A.  Table  33  shows  the  phase  speeds  for 

waveguides  of  5.3  cm  and  5.045  cm  (c_,.  and  c_„)  based  upon  a  free  field 

pr 

sound  speed  of  345  m/s.  Also  listed  is  the  expected  relative  phase 
shift  for  each  frequency  for  which  data  is  available.  Table  34  lists 
the  observed  relative  phase  shifts  (A<£obn),  the  average  (<4<£obn>),  and 
standard  deviation  for  each  frequency.  Tables  35  through  38  list  some 
possible  values  for  cpr,  the  phase  speed  in  the  rough  waveguide, 
calculated  for  cQ  =  345  m/s  and  A<^obn  +  360q  degrees  for  various  values 
of  q.  Also  listed  are  the  percent  change  in  epr.  It  is  seen  that  the 
randomly  rough  waveguide  seems  to  have  a  greater  effect  on  phase  speed 
than  did  the  deterministic  surface. 

An  additional  consideration  with  the  randomly  rough  surface 
in  the  waveguide  is  dispersion  due  to  attenuation.  The  attenuation  in 
the  randomly  rough  waveguide  is  considerably  more  than  that  in  the 
smooth  waveguide  at  the  lower  frequencies  of  7750  Hz  and  15750  Hz. 

Using  the  equations  of  section  V-A,  it  is  seen  that  dispersion  due  to 
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Table  32 


Pressure  Spectrum  Magnitude  (mV)  as  a  Function  of  Depth 
for  Raw  Data,  Mode  2  Driven,  Randomly  Rough  Surface. 


Signal-to-Noise  ratio  >  20dB  for  all  data. 


Data  at  31250  Hz. 


Depth 

1 .8 

Range 

2.7 

(meters) 

3.7 

4.6 

5.5 

0 

55.6 

12.01 

7.31 

4.52 

2.40 

1 

52.4 

12.89 

6.90 

4.65 

2.56 

2 

31.3 

14.12 

6.51 

4.26 

2.65 

3 

18.5 

7.85 

5.61 

3.42 

2.15 

4 

38.6 

9.96 

2.57 

3.18 

1 .99 

5 

53.4 

11.34 

2.58 

2.63 

1 .46 

6 

45.1 

10.93 

2.04 

1.63 

0.96 

7 

52.8 

11.30 

2.16 

0.49 

0.66 

8 

18.1 

4.39 

1  .96 

0.56 

0.65 

9 

43.0 

9.88 

0.73 

0.69 

0.56 

Data  at 

23500  Hz. 

Range 

(meters) 

Depth 

1 .8 

2.7 

3.7 

4.6 

5.5 

0 

188 

29.8 

10.08 

8.20 

2.06 

1 

148 

35.3 

10.10 

4.86 

2.22 

2 

110 

20.4 

6.14 

4.56 

2.46 

3 

29 

7.6 

4.09 

2.33 

2.03 

4 

106 

10.5 

4.98 

2.09 

1 .38 

5 

163 

22.7 

6.87 

2.81 

0.99 

6 

206 

36.0 

7.82 

2.96 

0.92 

7 

171 

28.5 

5.06 

1 .67 

0.76 

8 

39 

14.7 

1 .30 

0.34 

0.57 

9 

110 

15.6 

3.90 

1 .84 

0.82 
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Table  31  (eont'd) 

Mode  2  at  23500  Hz 


Smooth  Random  Roughness 


energy 

in  mode 

(X) 

energy 

in 

mode 

(X) 

Range 

(m) 

0 

1 

2 

oth 

0 

1 

2 

3 

oth 

1.8 

7 

4 

79 

10 

4 

27 

46 

11 

12 

2.7 

11 

6 

77 

6 

6 

43 

39 

6 

6 

3.7 

11 

19 

67 

3 

10 

31 

53 

2 

4 

4.6 

18 

4 

74 

4 

42 

38 

3 

9 

8 

5.5 

14 

6 

74 

6 

36 

54 

3 

0 

7 

Mode  2  at  31250  Hz 


Smooth  Random  Roughness 


energy 

in 

mode 

(X) 

energy 

in 

mode 

(X) 

Range 

(m) 

0 

1 

2 

3 

oth 

0 

1 

2 

3 

oth 

1.8 

11 

1 

72 

5 

11 

5 

36 

41 

7 

11 

2.7 

11 

8 

72 

2 

7 

5 

53 

30 

3 

9 

3.7 

3 

5 

68 

12 

12 

9 

74 

9 

2 

6 

4.6 

3 

6 

76 

8 

9 

34 

47 

8 

1 

10 

5.5 

8 

17 

66 

1 

4 

56 

41 

1 

0 

2 

Table  31 


Energy  Distribution  for  Mode  2  Driven  System, 
Smooth  and  Randomly  Rough  Surfaces. 

Modal  energy  distribution  by  percent  for  ten  depth 
measurements  taken  by  the  Nicolet  wave  analyzer 
and  processed  on  the  IBM  3030. 

Signal-to-noise  ratio  greater  than  20  dB  for  all  data. 

Mode  2  at  7750  Hz 


Smooth  Random  Roughness 

energy  in  mode  (X)  energy  in  mode  (X) 


Range 

(m) 

0 

1 

2 

oth 

0 

1 

2 

oth 

1  .8 

23 

1 

76 

0 

8 

3 

84 

71 

2.7 

24 

4 

72 

0 

12 

19 

75 

35 

3.7 

11 

3 

85 

1 

15 

7 

54 

24 

4.6 

22 

20 

57 

1 

86 

10 

3 

1 

5.5 

31 

13 

54 

2 

86 

9 

4 

1 

6.4 

40 

9 

51 

0 

93 

4 

3 

0 

Mode  2  at  15750  Hz 


Smooth  Random  Roughness 

energy  in  mode  (X)  energy  in  mode  (X) 


Range 

(m) 

0 

1 

2 

oth 

0 

1 

2 

oth 

1.8 

13 

2 

83 

2 

1 

3 

85 

11 

2.7 

11 

1 

84 

4 

4 

4 

82 

10 

3.7 

11 

6 

79 

4 

3 

1 1 

80 

6 

4.6 

12 

7 

80 

1 

14 

9 

61 

16 

5.5 

12 

3 

83 

2 

24 

51 

20 

5 

6.4 

16 

1 

81 

2 

38 

42 

17 

3 

Table  30  (cont’d) 


Mode  1  at  23500  Hz 


Smooth  Random  Roughness 


energy 

in 

mode 

(X) 

energy 

in 

mode 

1 

2 

3 

oth 

0 

1 

2 

3 

69 

8 

19 

2 

36 

7 

52 

3 

70 

1 

20 

5 

41 

50 

5 

2 

57 

2 

31 

3 

49 

30 

18 

0 

72 

1 

22 

4 

53 

40 

2 

2 

63 

0 

30 

3 

51 

43 

1 

1 

64 

2 

26 

7 

65 

31 

3 

1 

64 

5 

27 

2 

64 

25 

8 

1 

61 

2 

30 

3 

66 

25 

5 

0 

Mode  1  at  31250  Hz 


Smooth 

energy  in  mode  (X) 


Random  Roughness 
energy  in  mode  (X) 


Table  30 


Energy  Distribution  for  Mode  1  Driven  System, 

Smooth  and  Randomly  Rough  Surface. 

Modal  energy  distributions  by  percent. 

All  data  processed  by  Nicolet  using  ten  depth  measurements. 
Signal-to-noise  ratio  greater  than  20  dB  for  all  data. 


Mode  1  at  7750  Hz 


Smooth  Random  Roughness 

energy  in  mode  (.%)  energy  in  mode  (%) 

Range  0  1  2  oth  0  1  2  oth 

(m) 

1.8  0  96  0  4  1  89  9  1 

2.7  2  94  2  2  0  99  1  0 

3.7  1  97  2  0  1  98  1  0 

4.6  1  99  0  0  1  94  5  0 

5.5  1  99  0  0  5  92  2  1 

6.4  1  99  0  0  4  90  6  0 

7.3  0  100  0  0  2  89  6  3 

8.2  1  98  1  0  4  88  6  2 


Mode  1  at  15750  Hz 


Smooth  Random  Roughness 


energy  in 

mode 

(%) 

energy  in 

mode 

(X) 

Range 

(m) 

0 

1 

2 

oth 

0 

1 

2 

oth 

1.8 

1 

75 

3 

21 

12 

79 

1 

8 

2.7 

3 

73 

3 

21 

22 

75 

1 

2 

3.7 

2 

71 

0 

27 

24 

71 

1 

4 

4.6 

2 

74 

0 

24 

33 

58 

6 

3 

5.5 

1 

40 

0 

59 

15 

76 

5 

4 

6.4 

3 

84 

1 

12 

43 

47 

6 

4 

7.3 

2 

91 

1 

6 

20 

68 

5 

7 

8.2 

1 

98 

0 

1 

34 

53 

7 

6 
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The  attenuation  computed  using  the  equation  for  AI  is  frequency 
selective.  For  either  mode  1  or  2,  the  attenuation  rate  at  23500  Hz  is 
greater  than  the  rate  at  31250  Hz.  If  the  attenuation  were  a  function 
of  kh  for  kh<1,  Rayleigh  scatter,  attenuation  would  be  expected  to 
increase  with  frequency.  For  the  rms  height  of  1.5  mm,  kh  is  tabulated 
as  a  function  of  frequency  (Table  29)  and  is  shown  to  be  in  the  region 
of  Rayleigh  scatter  (kh<1): 

Table  29 

kh  As  a  Function  of  Frequency,  h=1.5  mm 
Frequency  (Hz) 

7750  15757  23500  31250 

kh  0.21  0.43  0.64  0.86 

The  coupling  of  energy  to  other  modes  is  almost  exclusively  to 
lower  order  modes  (Tables  30  and  31)  and  a  dramatic  oscillation  of 
percent  energy  among  modes  with  range  is  not  observed  with  the  randomly 
rough  surface.  Energy  is  almost  always  removed  from  the  driven  mode  at 
each  successive  range  except  for  the  mode  1  driven  system  at  15750  Hz 
and  the  mode  2  driven  system  at  23500  Hz.  The  reason  for  these 
exceptions  is  unclear.  The  mode  1  driven  system  couples  most  of  its 
energy  to  mode  0,  but  the  mode  2  driven  system  divides  the  energy 
between  modes  0  and  1  about  equally  for  frequencies  from  15750  Hz  and 
up.  An  interesting  point  about  the  mode  2  driven  system  can  be  seen  in 
the  raw  data  (Table  32).  As  the  range  from  the  source  increases,  the 
magnitude  of  the  observed  pressure  is  always  a  minimum  on  the  rough 
surface  at  the  bottom  of  the  waveguide.  This,  along  with  the  near  even 
distribution  of  percent  energy  in  modes  0  and  1,  suggests  that  the 
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Figure  22d.  Energy  Attenuation  at  31250  \lz.  Mode  1  Source, 
Smooth  vs.  Randomly  Rough  Surface. 
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Figure  22b.  Energy  Attenuation  at  15750  Hz,  Mode  1  Source 
Smooth  vs.  Randomly  Rough  Surface. 


112 


TABLE  23 


i 

Randomly  Rough  vs  Smooth  Surface  Attenuation  Rates,  Mode  1  Driven 


Mode  1  Energy  Attenuation 


Freq  (Hz) 

Smooth  (dB/m) 

Rough  (dB/ 

7500 

0.5 

3.2 

15750 

1.6 

4.8 

23500 

2.3 

5.1 

31250 

2.5 

4.4 

Mode  2  Energy  Attenuation 

Freq  (Hz) 

Smooth  (dB/m) 

Rough  (dB/ 

7750 

2.0 

4.8 

15750 

2.8 

10.5 

23500 

3.5 

12.2 

31250 

4.3 

10.1 

I 

I 


q 

cpr.  (m/s) 

_PE - Ell  x  100(2) 

c 

Pr 

-2 

340.68 

3.61 

-l 

348.93 

1.27 

0 

357.62 

-1.18 

+i 

366.72 

-3.76 

Table  36.  Change  in  Phase  Speed  of  Mode  1  over  Random 
Roughness  at  15150  Hz  Due  to  Cycle  Error  q. 
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q 

Cpr *  (m/ s ) 

Pr-~  Cpr'  x  100(%) 

C 

pr 

-3 

335.50 

3.79 

-2 

340.82 

2.26 

-1 

346.31 

0.69 

0 

352.00 

-0.94 

1 

357.86 

-2.63 

2 

363.92 

-4.36 

Table  37*  Change  in  Phase  Speed  of  Mode  1  over  Random 
Roughness  at  23500  Hz  Due  to  Cycle  Error  q. 


attenuation  at  these  frequencies  can  cause  a  change  in  phase  speed  of 
0.2**%  for  7750  Hz  and  O.H%  for  15750  Hz.  The  shift  is  toward  a  slower 
phase  speed  in  the  randomly  rough  walled  waveguide  due  to  the  increased 
viscous  wall  losses.  At  the  higher  two  frequencies,  the  difference  in 
attenuation  in  the  rough  versus  smooth  walled  waveguide  can  only 
account  for  a  slowing  of  the  phase  speed  in  the  rough  waveguide  of  less 


VI.  SUMMARY 


This  experiment  was  the  preliminary  investigation  of  mode 
interaction  due  to  steep  sloped  roughness  elements.  It  demonstrated 
the  possibility  of  obtaining  effective  cylindrical  spreading  in  a  long 
rectangular  waveguide  with  absorbing  sides.  Using  effective  single 
mode  transduction,  it  was  shown  that  attenuation  in  the  waveguide  was 
different  for  different  modes.  Measurements  of  the  interaction  of  low 
order  modes  due  to  the  existence  of  steep  sloped  roughness  demonstrated 
that  the  degree  of  mode  interaction  in  a  rough  walled  waveguide 
increases  with  increasing  Rayleigh  scattering  parameter  kh.  It  was 
also  demonstrated  that  the  existence  of  roughness  causes  a  significant 
change  in  the  phase  speed  of  the  modes  investigated.  The  precise 
amount  of  this  change  and  any  trends  with  frequency  have  yet  to  be 
determined.  With  mode  two  driven  over  the  randomly  rough  rigid 
surface,  the  data  indicated  a  range  dependent  change  of  the  rough 
walled  boundry  condition  from  rigid  toward  pressure  release.  This  too 
is  an  area  which  needs  further  investigation. 


APPENDIX  A 


MODE  SOURCE  TRANSDUCERS  AND  DRIVER  SYSTEM 

Effective  transmission  of  a  single  mode  in  a  waveguide  [Ref.  7]  is 
possible  if  the  driving  elements  are  located  to  reinforce  the  mode 
pressure  maxima  with  the  phase  between  adjacent  segments  shifted  130° 
(Fig.  A-1).  Using  this  principle,  an  amplifier  was  designed  to  split  a 
single  input  into  five  independent  channels  with  separate  amplitude  and 
phase  control  for  each  channel.  Phase  selection  between  channels  is 
either  0°  or  180°.  However,  at  low  amplifier  gain,  the  phase  between 
reversed  channels  is  only  170°,  as  observed  on  a  dual  trace 
oscilloscope.  This  is  the  maximum  phase  error  detected  and  is  not 
observed  at  high  amplifier  gain.  The  design  is  implemented  (Fig.  A-2) 
using  747C  dual  stage,  general  purpose  operational  amplifiers.  The 
first  stage  of  each  chip  is  a  buffer  to  isolate  that  channel  from  the 
others;  the  second  stage  provides  amplitude  control.  Phase  is 
controlled  by  a  toggle  switch  for  each  channel. 

Construction  of  the  mode  transducers  is  straightforward.  The 
transducer  measures  5  by  10  cm  which  covers  the  entire  depth  (5  cm)  of 
the  waveguide.  The  guiding  design  principle  is  to  have  no  element 
cross  a  null  of  the  desired  mode  and  to  maximize  the  element  area.  The 
larger  the  area,  the  stronger  the  element's  source  level.  For  a  mode 
two  source,  three  elements  are  imbedded  in  the  lucite  backing  so  that 
they  abut  at  the  two  nulls  for  mode  2,  1/4  depth  and  3/4  depth  (Fig. 
A-3).  The  aluminum  elements  are  bonded  to  the  lucite  backing  by  epoxy. 
The  unit  is  then  ground  to  make  the  face  of  uniform  thickness,  but 
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slightly  rough.  When  mylar  is  stretched  over  this  face,  the  roughness 
creates  an  air  space  between  the  face  and  the  mylar,  improving  the 
transducer  response  [Ref.  81.  The  mylar  is  pre-tensioned  in  a  jig  to 
apply  a  uniform  tension  across  the  mylar  sheet.  The  backplate  then 
receives  a  light  layer  of  quick  bonding  glue  about  its  perimeter  and  is 
pressed  against  the  mylar  sheet  in  the  jig.  After  30  seconds  the  glue 
has  bonded  and  the  mylar  is  released  from  the  jig.  A  ground  lead, 
common  with  the  center  tap  on  the  amplifier,  is  connected  to  its 
aluminized  side  and  provides  the  zero  plane  for  phase  reference 
(direction  of  mylar  displacement).  Curves  of  the  driving  response  for 
each  element  show  that  this  tensioning  technique  is  results  in  uniform 
response  (Fig.  A-4). 

The  transducer  requires  a  DC  bias  for  proper  operation.  This  bias 
is  provided  by  an  external  biasing  circuit  (Fig.  A-5)  which  can  be 
monitored  to  ensure  consistent  transducer  response  for  each  data  run. 
Initially,  attempts  were  made  to  permanantly  polarize  the  mylar  with  an 
electric  field  of  17,000  volts/cm.  The  polarization  uniformity  over  a 
large  area  of  mylar  was  uncertain,  however,  and  had  the  undesirable 
property  of  depolarizing  with  age  [Ref.  8].  Solving  this  stability 
problem  would  greatly  simplify  the  experimental  set-up  and  provide 


easy-to-use  sources. 


RESPONSE  FOR  EACH  ELEMENT 


2)  Wavetek  Phaselock  Signal  Generator  (WTK) 

a)  Mode  -  gated 

b)  Attn  -  dB 

c)  Sign  -  sinusoid 

d)  Freq  -  as  required 

C.  Polarization  Voltage  -  230  VDC 

D.  5  Channel  Amplifier  (Voltages  measured  on  RCA  WV-98c) 

1)  Mode  1  source  connections 

a.  Channel  1-7.2  volts  peak-to-peak,  switch  up 

b.  Channel  2-7.5  volts  peak-to-peak,  switch  down 

2)  Mode  2  source  connections 

a.  Channel  1-7.0  volts  peak-to-peak,  switch  up 

b.  Channel  2  -  4.0  volts  peak-to-peak,  switch  down 

c.  Channel  3  -  7.6  volts  peak-to-peak,  switch  up 

E.  PARS 

1 )  Gain  -  50 

2)  LF  roll  off  -  100  Hz 

3)  HF  roll  off  -  100  kHz 

F.  Kron-Kite  Filters 

1 )  High  pass  -  6  kHz 

2)  Low  pass  -  41  kHz 

G.  Nicolet  FFT  analyzer 

1 )  Signal  enhance 

2)  Average  Time 

3)  Transient 

4)  50-100  averages  -  AVG  N 

H.  Decade  Frequency  Synthesizer  -  249,999.9  Hz 
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APPENDIX  C 


EQUIPMENT  SETTINGS 

A.  Timing  Simulator 

1)  When  ARB  used  as  signal  generator 


Word 

Output 

Time 

0 

n.a. 

100  msec 

1 

1  (to 

ARB) 

100  usee 

2 

n.a. 

115  msec 

3 

9  (to 

OPHELEA) 

5  msec 

4 

n.a. 

100  msec 

2)  When  Wavetek  used  as  signal  generator 


Word 

Output 

Time 

0 

n.a. 

100  msec 

1 

1  (to 

WTK) 

1  msec 

2 

n.a. 

115  msec 

3 

9  (to 
& 

OPHELEA 

WTK) 

5  msec 

4 

n.a. 

100  msec 

B.  Signal  Generators 

1)  Abitrary  Waveform  Generator  (ARB) 

a)  Mode  -  trig 

b)  Output  -  on 

c)  amp  -  2.3  V 

d)  Pr3t  length  -  550 

e)  Block  -  full 

f)  Clock  -  internal 

g)  func 

i.  Composite  Signals  -  use  block  rate  of  3906  Hz 
RAM  8 

X  Y 

0  0 

16  122 

32  0 

ii.  Sinusoidal  Signals  -  use  desired  block  rate 
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Wavetek  Model  144  HF  Sweep 
Generator 


PMC  Regulated  Power  Supply 


Hewlett  Pcakard  Model  721A 
Power  Supply  (1) 


Hewlett  Packard  Model  467A 
Power  Amplifier  (1) 


Hewlett  Packard  Model  467A 
Power  Amplifier  (2) 


Power  Design  Inc.  (PDI) 
Model  3650R  Transistorized 
Power  Supply 


Hewlett  Packard  Model  721 A 
Power  Supply  (2) 


20  KHz  Square  Wave 


10  Volts,  Current  Limit  Fully 
Clock-wise 


5  Volts,  225mA  Short  Circuit 
Current 


XI 


24  Volts,  4  Amp  Current  Limit 


14  Volts,  225mA  Short  Circuit 
Current 


the  number  345  would  indicate  pin  5  of  the  chip  in  the  third  row  and 
fourth  column.  A  zero  indicates  that  the  pin  is  not  connected. 

Entries  at  pin  locations  which  are  not  numbers  refer  to  switches  indi¬ 
cated  in  Figure  B-5.  SQW  stands  for  the  output  of  the  20  KHz  square 
wave  generator. 

In  addition  to  the  above  mentioned  equipment,  four  commercial  DC 
voltage  supplies  are  used  to  run  the  device.  Two  power  amplifiers  are 
also  needed  to  step  up  the  output  current  of  the  timing  simulator  for 
triggering  the  stepping  motor  and  closing  the  solenoid  relay.  A 
diagram  of  the  equipment  setup  is  given  in  Figure  B-7.  Equipment 
identification  and  settings  are  listed  in  Table  B-1. 
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Figure  B-6.  Pin  Connects  for  Logic  Circuit. 
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possible  needle  positions  plainly  visible.  A  close  up  of  the  wire 
wound  resistor  and  sliding  contact  is  shown  in  Figure  B-3.  Figure  3-4 
shows  the  solenoid  which  raises  the  arm.  The  wire  connecting  the 
solenoid  to  the  lower  part  of  the  arm  is  visible. 

The  movement  of  the  arm  and  stepping  motor  and  data  sampling  by  the 
computer  are  sequentially  triggered  by  a  logic  circuit  (Fig.  B-5) 
driven  with  a  20  KHz  square  wave.  The  16  D  type  flip-flops  are  simply 
used  to  frequency  divide  the  20  KHz  square  wave.  (During  initial 
development  of  the  RSM  a  commercial  timing  simulator  was  used  in  place 
of  the  logic  circuit  and  square  wave  generator.)  It  should  be  noted 
that  this  timing  circuit  has  a  secondary  purpose  of  driving  a  device 
not  used  in  this  experiment.  Consequently  additional  logic  gates  and 
switches  are  evident  in  Figure  B-5  which  are  not  relevant  to  the 
operation  of  the  RSM. 

The  logic  circuit  was  constructed  on  standard  breadboards  with  the 
logic  chips  arranged  in  four  rows.  The  pin  connections  for  the  chips 
are  given  in  Figure  B-6.  Each  chip  is  assigned  a  two  digit  number 
which  describes  its  location  on  the  breadboard.  The  first  digit 
represents  the  row  the  chip  is  in,  the  second  digit  represents  its 
column.  This  number  is  located  in  the  center  of  each  chip  in  the 
Figure.  Below  this  number  is  a  number  indicating  the  chip  type.  Each 
chip  has  14  pins.  Pin  1  is  in  the  lower  left  corner  for  each  chip. 

Pins  are  numbered  successively  counter  clock-wise  around  the  chip.  The 
numbers  at  pin  locations  indicate  to  which  pin(s)  that  pin  is 
connected.  The  first  two  digits  indicate  the  chip,  the  third  (and 
fourth)  digit(s)  indicate  the  pin  number  of  that  chip.  For  example. 


for  placement  of  the  needle  for  measuring  different  ranges  of  rough¬ 
ness.  The  low  end  of  the  range  is  determined  by  the  minimum  height 
deviation  detectable.  The  high  end  of  the  range  is  determined  by  the 
maximum  displacement  of  the  needle  which  will  ensure  that  the  sliding 
contact  does  not  lose  contact  with  the  resistor.  This  is  set  by  the 
ratio  of  length  of  the  long  arm  to  length  from  the  pivot  point  to  the 
needle  position.  With  the  needle  in  the  position  to  provide  the 
greatest  resolution,  roughnesses  in  the  range  of  0.02mm  to  2.85mm  can 
be  measured.  With  the  needle  in  the  other  two  positions,  the  ranges  of 
roughness  which  can  be  measured  are  0.03mm  to  4.28mm  and  0.05mm  to 
5.70mm.  Calibration  runs  were  made  on  a  flat  glass  plate  to  determine 
the  variance  of  the  readings.  Calibration  runs  along  a  smooth  aluminum 
wedge  were  made  to  determine  the  voltage  to  height  conversion  factors. 

After  a  point  is  sampled,  a  relay  is  closed  which  energizes  a 
solenoid  above  the  lower  section  of  the  arm,  raising  the  arm  by  a  steel 
wire.  A  stepping  motor  is  then  triggered  to  turn  a  threaded 
positioning  rod.  The  rod  moves  the  arm/needle/resistor  assembly  over 
the  next  point  to  be  sampled.  Horizontal  spacing  between  sampled 
points  is  0.0355mm.  With  only  minor  changes  in  the  control  circuitry, 
this  can  be  changed  to  any  integral  multiple  of  0.0071mm. 

The  surface  to  be  sampled  is  placed  on  a  flat  test  bed  made  of 
3/16  inch  thick  aluminum  plate.  This  plate  can  be  raised  or  lowered  by 
means  of  height  adjustment  screws  located  at  each  corner.  Downward 
deflection  of  the  test  bed  is  minimized  by  inflating  the  air  bladders 
beneath  it. 

Figure  B-2  is  a  close  up  of  the  needle  and  lower  arm  with  the  three 
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APPENDIX  B 


RANDOM  SURFACE  MEASURING  DEVICE 

The  random  surface  measuring  device  (RSM)  measures  the  height 
deviation  from  a  mean  plane  of  a  surface  (roughness)  for  roughnesses 
in  the  range  of  0.02mm  to  5.7mm  with  a  possible  accuracy  of  0.02mm. 
Height  measurements  are  taken  at  horizontal  intervals  of  0.0355mm. 

These  height  measurements  are  converted  to  proportional  voltages  which 
are  stored  directly  on  OPHELEA's  magnetic  disk  memory.  1024  points, 
the  maximum  possible  with  the  software  used,  are  recorded  on  each  data 
run  and  processed  to  give  statistics  of  the  surface  such  as  RMS  height, 
RMS  slope,  and  correlation  length.  The  statistics  of  several  runs  are 
then  averaged  to  produce  the  statistics  given  in  the  Section  IV  which 
describes  the  randomly  rough  surface. 

The  construction  of  the  RSM  is  shown  in  Figure  B-1.  A  constant 
voltage  is  applied  across  the  wire  wound  resistor  near  the  top  and  can 
be  adjusted  with  the  potentiometer  at  the  base.  The  "L"  shaped  arm  is 
the  key  component  with  a  needle  attached  under  its  short  section.  The 
long  section  has  an  electric  contact  at  its  end  that  slides  across  the 
wire  wound  resistor.  The  arm  pivots  freely  at  the  elbow  and  will 
rotate  by  gravity  until  the  needle  rests  on  the  rough  surface.  The 
voltage  difference  between  the  sliding  contact  and  the  positive  end  of 
the  resistor  is  proportional  to  the  distance  the  arm  has  moved  for  the 
needle  to  contact  the  surface.  The  voltage  is  thus  a  direct  measure  of 
the  height  at  that  point  of  the  surface. 


On  the  lower  part  of  the  arm,  there  are  three  possible  positions 
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